Dissipative vortex solitons in 2D-lattices 
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We report the existence of stable symmetric vortex-type solutions for two-dimensional nonlinear 
discrete dissipative systems governed by a cubic-quintic complex Ginzburg-Landau equation. We 
construct a whole family of vortex solitons with a topological charge S = 1. Surprisingly, the 
dynamical evolution of unstable solutions of this family does not alter significantly their profile, 
instead their phase distribution completely changes. They transform into two- charges swirl-vortex 
solitons. We dynamically excite this novel structure showing its experimental feasibility. 
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The study of discrete nonlinear systems has been an 
active area of research during the last twenty years due 
to its broad impact in diverse branches of science and to 
its potential for technological applications [H-Q. Until 
now, nonlinear optics has been the main scenario chosen 
to test this phenomenon, essentially because of both, its 
comparative experimental simplicity and its direct con- 
nection with theoretical models. Nonlinear self-localized 
structures, usually termed as discrete solitons, have been 
predicted and observed for one- and two-dimensional ar- 
rays [H|, Q . A discrete vortex soliton is denned as a non- 
linear self-localized structure whose phase changes 2irS 
radians azimuthally. S is an integer number known as 
the vorticity or topological charge of the solution. The 
existence of discrete vortex solitons in conservative sys- 
tems have been reported on several works [7]. For the 
continuous case, dissipative vortex soliton families have 
been found to be stable for a wide interval of S- values [i| . 
Very recently, symmetric stable vortices have also been 
predicted in continuous dissipative systems with a peri- 
odic linear modulation Q. 

Nowadays, dissipative models offer a more complete 
and realistic description of different physical systems. 
In conservative models, gain and loss are completely 
neglected and the dynamical equilibrium is reached by 
means of a balance between nonlinear and dispersive ef- 
fects. For dissipative systems, it must also exist an addi- 
tional balance between gain and losses, turning the equi- 
librium into a more complex process MOL The Ginzburg- 
Landau equation is - somehow - a universal model where 
dissipative solitons are their most interesting solutions. 
This model appears in many branches of science like, 
for example, nonlinear optics, Bose-Einstein condensates, 
chemical reactions, superconductivity and many oth- 
ers PH. 

In this work we deal with discrete vortex solitons in 
dissipative 2D lattices governed by a discrete version of 
the Ginzburg-Landau equation. We have found differ- 
ent families of these localized solutions connected succes- 
sively by means of saddle- node bifurcations. We studied 



their stability and found two types of stable vortex fam- 
ilies coexisting for the same set of parameters. We have 
dynamically unveiled the second type of stable solution 
by following the decaying of an initially unstable vortex. 
This observation is very different to the results shown 
in Ref. where an unstable vortex just vanishes on 
propagation, through completely radiative decay. More- 
over, our final vortex solution possesses a nontrivial phase 
structure where two different charges coexist. 

Beam propagation in 2D dissipative waveguide lattices 
can be modeled by the following equation: 

ilpm,n H~ Clpm,n + 1 1pm,n \^1pm,n H~ ^\lpm,n\^1pm,n ~ 
^1pm,n + IS 

\lpm,n\ 1pm, n (i) 

Eq.(pQ) represents a physical model for open systems that 
exchange energy with external sources and it is called 
(2 + 1) discrete complex cubic-quintic Ginzburg-Landau 
equation. ipm,n is the complex field amplitude at the 
(m, n) lattice site and ip m ,n corresponds to its first deriva- 
tive respect to the propagation coordinate z. The set 
{m = 1 , . . . , M} x {n = 1 , . . . , N} defines the array, being N 
and M the number of sites in the horizontal and vertical 
directions (in all our computations N = M = 17). The 
fields propagating in each waveguide interact only with 
nearest-neighbors through their evanescent tails. This 
interaction is described by the discrete diffraction oper- 
ator Clpm^n = C(VWl,n + 1pm — l,n + 1pm,n+l + 1pm,n — l)i 

where C is a complex number. Its real part indicates 
the strength of the coupling between different sites and 
its imaginary part denotes the gain or loss originated by 
this coupling. The nonlinear higher order Kerr term is 
represented by v while e > and \i < are the co- 
efficients for cubic gain and quint ic losses, respectively. 
Linear losses are determined by negative S. 

Unlike the conservative discrete nonlinear Schrodinger 
(DNLS) equation, the power defined as 
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Q(z)= IVwWI 2 
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is not a conserved quantity in the present model. How- 
ever, for a self-localized solution, the power and its evolu- 
tion will be the main magnitude that we will monitor in 
order to identify different families of stationary solutions. 

We look for stationary solutions of Eq.(pQ) of the form 
^m,n(z) = m , n exp[iAz] where m , n are complex num- 
bers and A is real; also we are interested that the phase 
of solutions change azimuthally an integer number (S) 
of 2tt. In such a case the self-localized solution is called 
a discrete vortex soliton [12| with vorticity S. By in- 
serting the previous ansatz into model (pp) we obtain the 
following set of algebraic coupled equations: 

— A0 m ,n + C0 m , n + |0ra,n| 2 0ra,n + ^|0ra,n | 4 0ra,n = 

iS</> m ,n + ^#m,n| 2 0m,n + ^|0m,n| 4 0m,n • (3) 

We look for vortex-type solutions by solving equations 
(|3]) with a multi-dimensional Newton-Raphson iterative 
algorithm. The method requires an initial guess that we 
construct as explained below. 

In the high- confinement limit, single peak solutions 
(fundamental bright solitons) were predicted to exist 
in dissipative nonlinear media [13|- In that limit, 
we obtain the following approximation: 0q w — (e + 
V / ? 7 ¥)/(2/i), A * 0§ + i/0g, and a w |C0 o /(A + ^)|. 
Here, 0o corresponds to the central amplitude, A the non- 
linear propagation constant, and a the first adjacent am- 
plitudes. We can see from the above that the amplitude 
of each peak is a function of £, \i and 5; if we set the last 
two parameters, the amplitude takes a ^-quadratic form 
with e as the bifurcation parameter. Now, we place a 
single peak approximation at each corner of a square sub- 
lattice L as a superposition of four fundamental bright 
solitons [14j : 
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where a = 2a. Now, we define a phase operator as 
®m,n = exp{i[arctan(— n/m)]}, and write our initial S = 
1 ansatz as 

0m, n -^m,n ' ^m,n • 

With this initial guess, we construct a family of 4-peaks 
symmetric vortex solitons with vorticity 5 = 1, whose 
stability is monitored through a standard linear stabil- 
ity analysis [15] . Figure QJ a) shows a Q versus e dia- 
gram for these solutions including their stability. This 
figure shows the coexistence, for the same set of param- 
eters, of two different branches of stable solutions and, 
also, three different families of unstable solitons. Dif- 
ferent families are successively connected by saddle-node 
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Figure 1. (Color online) (a) Q versus s diagram for discrete 
vortex solitons. Continuous and dashed lines correspond to 
stable and unstable solutions, respectively, (b) and (c) Color 
map plots for the amplitude and phase profiles, respectively, 
of the solution for e = 1.1 indicated by a black dot in (a). 
C = 0.8, 5 = -0.9, a = -0.1, v = 0.1. 

bifurcation points. An example for a solution of branch 
A is shown in FigsOJb) an d (c) [black dot in FigQJa)]. 

This solution is very similar to our initial ansatz 
sketched in Eq.(jU) with a full topological charge S = 1. 
This agreement validates the seed we constructed as a 
first approach to find stationary vortex- type solutions. 
As the nonlinear amplification is diminished, the stable 
branch A reaches a first saddle-node point for e « 0.637. 
At this point, this family turns around and a new fam- 
ily emerges: the unstable branch labeled B. After that, 
two more saddle-node points appear connecting the new 
branches B with C and, then, C with D. The unstable 
branch D is mostly hidden because it is located at the 
same region that the stable branch labeled E. Branches 
A — D preserve the vorticity S = 1 while the ampli- 
tude profiles change adiabatically. It is worth mention- 
ing that branches A, D and E also exist for higher e 
values, with the power increasing monotonically, as the 
high- confinement limit predicts. 

As said before, in FigJTJa), curves D and E are in- 
distinguishable. In order to see their differences more 
clearly, we plot a zoom in FigJSfa) of region Q ~ 81.2 for 
a narrow region around the gray arrow in the FigQJa). 
The first solution on branch E [black dot in Figj2fa)] 
was obtained dynamically; i.e., we numerically integrated 
Eq.(pQ) by using an unstable solution [gray dot in Fig|2fa)] 
as initial condition. Contrary to previous observations, 
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Figure 2. (Color online) (a) Zoom of branches D and E for a 
narrow region around the gray arrow in the Fig[TJa). (b) Nu- 
merical simulation of model (pQ) showing the power transition 
sketched by an cyan arrow in (a) for s — 1.1. 



for the evolution of unstable vortex solitons [9|, we no- 
ticed that the power Q makes one oscillation and then 
stabilizes very rapidly around a new equilibrium value 
[see FigJSJb)]. This new value was indeed very close to 
the initial one, but now it corresponds to a new station- 
ary solution that propagates stably by keeping the same 
amplitude profile but a different phase structure. We 
took this new solution as an initial guess in our Newton- 
Raphson scheme and we constructed the whole stable 
branch E shown in FigOJa). 

The amplitude profile for solutions corresponding to 
the gray and black points at Fig|2fa) is shown in Figj3^a). 
This profile is almost identical for both solutions and it 
corresponds to a new structure that we define as ''''swirl- 
vortex soliton" . However, both solutions have a quite 
different phase profile. The unstable solution (belonging 
to branch D) possesses a full phase profile with charge 
5 = 1 [see Figl2fb)]. A very interesting thing related 
with charges happens with the stable swirl- vortex soliton 
[see Fig|3fc)]. For the first square contour [the innermost 
discrete square trajectory on the plane (n, m)] we can see 
that the vorticity has a 5 = 1 value, while for the next 
contours the vorticity has decreased to 5 = —3 [5 > 
(5 < 0) means a clockwise phase structure from — ir to 
7r (tt to — 7r)]. Therefore, there is a stable coexistence of 
two different topological charges for the same mode. This 
new type of structure would correspond to a " two- charges 
swirl-vortex soliton" and - as far as we know - this would 
be the first time they are predicted in nonlinear lattices. 

In order to go deeper in the understanding of this stabi- 
lization process and its dependence with the phase struc- 
ture, we found another example in which the change in 
vorticity is also related to the stabilization of the solution. 
As an initial ansatz, we constructed (in the same way as 
the four-peaks vortex) a symmetrically-centered twenty 
peaks 5 = 1 configuration. From the Newton-Raphson 
scheme we obtain an unstable vortex solution with the 
amplitude and phase profiles shown in Figs|4]^a) and (b), 
respectively. Again, we use this solution as an initial 
condition and numerically integrate Eq. (pQ). As a conse- 



(a) 


O 


O O 


O 


o 


o o 


o 




4-2 


2 4 



(b) 






t 





Figure 3. (Color online) Color map plots for the solutions 
indicated with dots in Fig|2ja). (a) Amplitude profile, (b) 
and (c) Phase profiles for the gray and the black points, re- 
spectively. 



quence of the larger number of excited sites, the power 
of this solution is higher, namely Q = 202.43. Similar to 
the swirl- vortex case, the solution evolves and converges 
to a stable solution with a very similar amplitude profile, 
but with a very different phase structure. For this case, 
the topological charge of the solution transforms from 
5 = 1 into 5 = — 3. It is worth to mention that there 
is a mismatch between the vorticity of the first and next 
contours of the lattice [see FigH^c)]. 
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Figure 4. (Color online) Color map plots for the twenty peaks 
vortex solution, (a) Amplitude profile, (b) and (c) Phase pro- 
files for the unstable 5=1 and the stable 5 = — 3 solutions, 
respectively. C = 0.8, 5 = —0.9, \i = —0.1, v = 0.1 and 
e = 1.1. 

Looking at the colormaps for the stable vortex soliton 
shown in FigsOJb) an d (c) we can realize that amplitude 
and phase structures have the same reflection and rota- 
tion symmetries. Similar happens with Figs |3] and FigsJU 
As it was shown in previous works the stability for one 
solution with a high number of excited sites requires an 
increment of its topological charge [HI, [l6[ . From this we 
understand why the dynamic evolution modify the vor- 
ticity of our solutions. So, we may conclude that the in- 
stability for complex-structure solutions of charge 5 = 1 
is essentially related with the geometric distribution and 
the number of excited sites. 

We have also explored the above phenomenology in 
conservative-cubic systems (DNLS limit): 5 = /jl = v — 
e = 0. There, we found two branches of swirl- vortex soli- 
tons; one with charge 5 = 1 and one with "two-charges" . 
The first one is always unstable while the two-charges 
solution is only stable for higher level of power. Here, 
we did not observe any decaying mechanism essentially 
because the system has not gain to increase the power 
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Figure 5. (a), (b), (c) and (d) Phase (upper) and amplitude 
(lower) profiles for different values of z marked by red points 
in (e) where the power evolution is displayed. 

and stabilize the solution. 

Finally, and with the aim of proposing a possible ex- 
perimental realization, we numerically integrate model 
(PQ) by taking - as an initial condition - a profile with 
eight peaks spatially distributed in the form of a swirl- 
vortex, including its phase structure S = 1, as FigEfa) 
shows. With this configuration we initialize the dynamic 
evolution observing that the system rapidly converges to 
an stable stationary two- charges swirl-vortex soliton [see 
Figs|5fa)-(d)]. We can see that the amplitude profiles 
slightly change during the propagation: the eight initial 
peaks remain almost unaltered. On the other hand, from 
FigsJSJd) we see that for the first square contour the vor- 
ticity is preserved being S = 1, while for the next con- 
tours the charge has transformed into S = — 3. FigsJSfe) 
shows the evolution of power with some initially small 
oscillations and, lately, a tendency to the stabilization of 
the profile. This example shows the robustness of our 
prediction and its chances to be observed in real dis- 
sipative systems because the initial condition could, in 
principle, be easily implemented in current experimental 
setups. Another very interesting point is that the system 
naturally evolves to a "2-charges" structure. Our initial 
condition has an unstable phase structure which guaran- 



tees the decaying to another type of mode, but not nec- 
essarily to the one we are interested in; it could perfectly 
just be destroyed by the internal dynamics [9]. However, 
the system favors the excitation of a swirl- vortex solution 
which propagates stably for long propagation distances. 

In conclusion, our results reveal the existence of dis- 
crete vortex solitons in dissipative 2D-lattices. We have 
found stable and unstable vortices by performing differ- 
ent continuation methods. In particular we concentrated 
the study on a new type of stable structure, the so-called 
two- charges swirl-vortex soliton. We were able to dy- 
namically excite it by using a simple initial configuration 
and therefore, believe in the feasibility of experimental 
observation of this novel type of dissipative structures. 
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